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Abstract 



The "basketball diagram" is a three-loop vacuum diagram for a scalar field 
theory that cannot be expressed in terms of one-loop diagrams. We calculate 
this diagram for a massive scalar field at nonzero temperature, reducing it to 
expressions involving three-dimensional integrals that can be easily evaluated 
numerically. We use this result to calculate the free energy for a massive 
scalar field with a 4 interaction to three-loop order. 
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I. INTRODUCTION 



One of the obstacles to making progress in thermal field theory is that the technology 
for explicit perturbative calculations is underdeveloped. The formalism of thermal field 
theory is sufficiently complicated that there are often theoretical issues that are difficult to 
resolve without explicit calculations. An example is the gluon damping rate, which in the 
conventional perturbative expansion is plagued by problems involving gauge invariance. A 
formal solution to the problem by the resummation of hard thermal loops was presented by 
Braaten and Pisarski in 1990 [|TJ. However, the solution was not widely accepted until the 
leading order expression for the damping rate was calculated explicitly 0. 

Around 1994, there was a significant step forward in the calculational technology for 
massless theories. The first perturbative calculation in thermal field theory that was carried 
out to high enough order that the running of the coupling constant came into play was 
a calculation of the free energy of a massless scalar field theory to order g 4 by Frenkel, 
Saa, and Taylor in 1992 ||. In 1994, there were several other calculations of the free 
energy to fourth order in the coupling constant: a calculation by Coriano and Parwani 
for QED [ffl and completely analytic calculations by Arnold and Zhai [El] for a massless 
scalar theory (correcting an error in Ref. 0) and for QCD. Arnold and Zhai made a 
particularly significant contribution by showing how three-loop vacuum diagrams, such as the 
so-called "basketball diagram" labeled 2b in Fig. [l|, could be evaluated analytically. These 
analytic calculations were then quickly extended to order g 5 for massless scalar theories 
H|7|], abelian gauge theories [§|§, and nonabelian gauge theories |10|,|lTJ. These explicit 



calculations revealed that the weak coupling expansion has convergence problems whose 
severity had not previously been appreciated. 

Calculations in theories that include massive particles at nonzero temperature are more 
difficult, because there is a second scale in the problem. The calculational technology for 
massive theories is much less well-developed than that for massless theories. In addition 
to the obvious applications to theories with massive particles, calculations with massive 
propagators may also be useful for massless theories. In such theories, some of the most 
important thermal corrections have the effect of generating masses for the massless parti- 
cles. These corrections may be responsible for the poor convergence properties of the weak 
coupling expansion. One of the most promising methods for resumming these corrections 
in scalar field theories is "screened perturbation theory" proposed by Karsch, Patkos, and 



Petreczky ||12|| . This method involves adding and subtracting a mass term from the La- 
grangian and treating the subtracted term as a perturbation. The integrals encountered in 
the screened perturbative expansion are those of the corresponding massive theory. Screened 
perturbation theory has been applied to the free energy at the two-loop level, and it seems 
to dramatically improve the convergence of the perturbative series ||12|| . To determine how 
effective this method is in resumming the large perturbative corrections, it is essential to 
calculate higher order corrections explicitly. 

In this paper, we take a step forward in the calculational technology for massive field 
theories by evaluating the basketball diagram. The diagram cannot be evaluated analytically, 
but we reduce it to expressions that involve integrals that are at most three-dimensional and 
can easily be evaluated numerically. Using our result for this diagram, we calculate the free 
energy for the massive 4 field theory through next-to-next-to-leading order in the coupling 
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constant g 2 . 



II. BASKETBALL DIAGRAM 

The basketball diagram is the diagram labeled 2b in Fig. [l|. At nonzero temperature T, 
it involves a three-fold sum-integral over Euclidean momenta: 



Jba11 jpqr P 2 + m 2 Q 2 + m 2 R 2 + m 2 (P + Q + R) 2 + m 2 ^ 

The Euclidean four-momentum is P = (p,p4 = 27mT), where n is an integer, and its square 
is P 2 = p 2 + p\. The sum-integral f P represents the sum over the Euclidean energies and 
the integral over the spatial momentum: 

i = T ^j ■ (2) 

JP P4, J P 

We use dimensional regularization of the integral over p to cut off the ultraviolet diver- 
gences in the sum-integral. Our convention for the measure in the integral is 

{ 4tt J J (2tt) 3 ^ ' 1 ' 

where 3 — 2e is the number of spatial dimensions, 7 is Euler's constant, and fi is an arbi- 
trary momentum scale. The factor of fi 2e gives the dimensionally-regularized integral ([I]) 
dimensions of (energy) 4 . 



A. Decomposition into integrals 

In order to evaluate the sum-integral ([[]), we first reduce it to integrals that contain 
factors of the Bose-Einstein distribution function n(E) = \j[eP E — 1) with positive energy 
E. We use the method of Bugrij and Shadura [(0| which expresses the coefficients of the Bose- 
Einstein factors in terms of S-matrix elements at zero temperature. This strategy was also 
used by Frenkel, Taylor and Saa [Q to evaluate the massless basketball diagram. Although 
the derivation by Bugrij and Shadura is lengthy, their final result can be obtained by making 
some simple substitutions in the expression (fj). The sum- integrals over Euclidean momenta 
P = (p,Pi) are replaced by integrals over Minkowski momenta p = (po,p): f P — > —i J p - 
The Euclidean propagators are replaced by the Minkowski propagators of the real-time 
formalism: 

-m]— 2 — i (~ % -r-r- + n(\p \)27r5(p 2 - m 2 )) , (4) 

P z + m z \p z — m z + is J 

where p 2 = p\ — p 2 . The sum- integral ([Tp is then given by the real part of the resulting 
expression, which now involves a three-fold integral over Minkowski momenta. For some of 
the integration momenta, the energy appears in the argument of a Bose-Einstein distribution. 
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It is convenient to Wick-rotate the remaining integration momenta back to Euclidean space: 
J — > i J P , where the measure of the dimensionally regularized integral over the Euclidean 
momentum P is 



47T 



4-2e 



(5) 



Having carried out this procedure, the basketball diagram can be expressed as the sum of 
terms with different numbers of Bose-Einstein factors: 



'ball 



X(0) +4 J(1) + 6T (2) +4:Z: (3) 



(6) 



The term with four Bose-Einstein factors is purely imaginary and it vanishes when we take 
the real part. The other terms are 



J(0) 

j(D 



1 



1 



1 



pqr P 2 + m 2 Q 2 + m 2 R 2 + m 2 (P + Q + R) 2 + m 2 
1 1 1 



n5{p) 



qr Q 2 + m 2 R 2 + m 2 (P + Q + R) 2 + m 2 

1 



p2._ 



2«=B.U)U)/^ 



+ m 2 (P + Q + R) 2 + m? 



{P+Q) 2 = -{p+q) 2 -i£ 



= Re / n6(p) f nS{q) f nS{r) 



(-1) 



(p + q + r) 2 — m 2 + is ' 
where we have used the shorthand nS(p) = n(\po\)2n5(p 2 — m 2 ). 



(7) 
(8) 
(9) 
(10) 



B. Zero thermal factors 

We first consider the term , which is the integral for the zero-temperature basketball 
diagram. The poles in e for this diagram have been calculated by Kastening WM and by 
Chung and Chung [15|]. The finite terms in the diagram can be obtained by following the 
strategy used in Appendix B of Ref. to calculate the zero-temperature basketball diagram 
in three dimensions. 

The diagram is first Fourier-transformed to coordinate space, which reduces it to an 
integral over a single coordinate R: 

J(°) = / V\R) , (11) 
Jr 

where the potential V(R) is 

and K u (z) is a modified Bessel function. The measure for the integration over R is 



4 



d^R. 



ir \ Air J 

After integrating over angles in 4 — 2e dimensions, the integral reduces to 



(13) 



J(0) 



m 



.3c 
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(47T) 



IY2 



/■oo 

/ tft r 1+2e Kt_ e (2t) 

Jo 



(14) 



The £ — > region of the integral gives poles in e. The small-£ behavior of the Bessel function 
is given by the power-series expansion 



#i- e (2t) 



r(l - e)£ 



-l+e 



rn-e)r( 



■E 

j=o 



£ 2 i+ l- 



(j' + l)!r(j + l + e) j!r(j + 2-e) 



(15) 



In the integral over £ in (0), the poles in e come from the £ 5+6e , £ 3+6e , £ 3+4e , £ 1+66 5 
i _1+4e , and £~ 1+2e terms. We can calculate the poles analytically by multiplying each of 
these terms by an appropriate convergence factor and integrating over £. After these terms, 
with their convergence factors, are subtracted from the original integrand, the remaining 
integral is convergent for e = and can be evaluated numerically. We choose convergence 
factors of the form (e 8 *) n e~ 8< , where {e x ) n is the truncated power series for the exponential 
function: (e x ) n = YJI=qX 1 This convergence factor behaves like 1 + 0{t n+l ) at small £, 
and has the same exponential falloff as Kf(2t) at large £. The resulting expression for the 
integral is in the limit e — > 



/ ci£ £- 1+2e K 1 4 _i2£) 
Jo 



r 4 n 
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dt £- 1+2e e" 8 * 



-4+4e 



4£ 



-2+2e 



,2c 



r(i + e) 
r(2 - C ) 



2t 



2c 



e l + e 



; 2 e _ 2r(2 



r(3 



>_r(i + e) 



r(2 



fOO 1 f 

+ y o dt-lK}(2t)- 



-St 



16 



-(e M ) 4 + -(21og£ + 2 7 -l)(e M ) 



£ 4 



£ 2 



+2(2 log£ + 2 7 - |) + 6(2 log£ + 2 7 - l) 1 



(16) 



The first integral in ( |I~6D can be evaluated analytically in terms of gamma functions, and it 
reduces in the limit e — > to 



, 1 17 59- 36vr 2 821 
r(l + 6e)<!- + — + — — + 381 log 2 



6e 2 



12e 



24 

617T 2 



-90 log 2 2 - 216 log 3 2 12tt 2 log 2 - 27^(1) 



(17) 



where ip(z) is the digamma function. The numerical value of the second integral in ( |T6"D 
is 0.36106. Inserting these results into (0) and keeping terms through order e°, our final 
result is 
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j(o) 1 (V- 



2 



2 23 35 + tt 
(4vr) 6 Vm7 p + 3? + ~~ 2l " 

where the numerical value of the constant is Cq = 39.429. 



m 4 , (18) 



C. One thermal factor 

The expression (H) for can be written as 

J« = J sun (-m 2 ) / n(\p \)2n5(p 2 -m 2 ), (19) 

where I sun (P 2 ) is the integral for the "setting sun diagram" in the boson self-energy at zero 
temperature: 

/su ^ P ^ = isfi Q 2 + m 2 R 2 + m 2 (P + Q + R) 2 + m 2 ' ^ 
The integral over p in ([TJ]) can be written as 

/ /(W )2^- m S) = ^(£) 2 V*, (21) 

where Ji is a function of /5m defined by ( |A.5| ) in the Appendix. 

The setting-sun integral (|20|) at P 2 = —m 2 can be evaluated by following the strategy 
used in Appendix B of Ref. to calculate the corresponding integral in three dimensions. 
After Fourier-transforming, it reduces to an integral over a single coordinate R: 



Isun ( TTl 



2', 



/ e iPR V 3 (R) 
Jr 



(22) 

P 2 =— m 2 



After averaging over angles in 4— 2e dimensions and evaluating at P 2 = — m 2 , the exponential 
factor becomes 



P 2 =—m 2 



r(2 " e) (^?) 1 e/l ^ (mjR) ' (23) 



where I v (z) is a modified Bessel function. The integral over R then reduces to a one- 
dimensional integral, and ( p2]) becomes 



9s m 2 / u x " 
-wr) ' 



(4tt) 



' 16e 27e dt t- 1+2e h_ e (2t)Kt e (2t). (24) 



The t — > region of the integral gives poles in e. The small-t behavior of the Bessel function 
Ii- e (2t) is given by the power-series expansion 

oo j.2j + l — e 
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In the integral over t in (|24]), the poles in e come from the t~ 3+4e , t~ 1+4e , and t~ 1+2e terms. 
We can calculate the poles analytically by multiplying each of these terms by an appropriate 
convergence factor and integrating over t. After these terms, with their convergence factors, 
are subtracted from the original integrand, the remaining integral is convergent for e = 
and can be evaluated numerically. We choose convergence factors of the form (e 6 *) n e~ 6< . 
The resulting expression for the integral is in the limit e —>■ 



/ dt t~ 1+2e h- e (2t) Kf(2t) 
Jo 



T 2 (l 



dt r 1+2e e- 6t 



-2+2e 



(«"),+ 



t 



2c 



+ - It* - 



f°° 1 f 

+ / dt -\l x (2t)Kl(2t) - 
J o t 



2-e e 

-6t r 1 1 
J- / J- 



£(1+6) 

r(2-c) 



/2V e 6 *) 2 + | + 3(2 logt + 2 7 -l) 



(26) 



The first integral in (p6|) can be evaluated analytically in terms of gamma functions, and it 
reduces in the limit e — > to 



rCl + ^l-A--- — + -1 OS6 + 31 og 2 6 + ^ 
V ; | 4e 2 8e 16 2 5 5 2 



(27) 



The numerical value of the second integral in ( p26|) is —1.2713. Inserting these results into 
(E£jj) and keeping terms through order e°, we obtain 



1 



-m 



(4vr) 4 U/ 



3 17 
"2^"47 + Cl 



(28) 



where the numerical value of the constant is C\ = —9.8424. Inserting (j2~8D and ( 21f) into 
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(|i9|), our final result is 

j(D 



(47T) 



6c 



17 + r 



2rTi2 



J im z T 



(29) 



Note that the integral J\ depends on e. 



D. Two thermal factors 

The expression (^) for X^ 2 ) involves the "bubble integral" 

W(P2) = S R R 2 lm 2 (P + R) 2 + m 2 > (3 ° } 
which can be evaluated using a Feynman parameter: 



1 ( I 1 \ 2e [l f 1 

^bubble (-P 2 ) = TT~v?\~) / dxlog 

(Ati) 1 \mj e Jo 



m 2 + x(l — x)P 2 
m 2 



(31) 



The real part of the integral evaluated at P 2 = —(p + q) 2 —ie is obtained by simply replacing 
the argument of the logarithm by its absolute value. When ( |3~T1) is inserted into @, the 
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coefficient of 1/e can be evaluated using (pT|). Reducing the other term to an integral over 
spatial momenta, we obtain 



1(2) 



(47r) 6 \m) 



^y e \-J 2 1 T 4 -32 I dp 



p 2 n(E p ) 
E n 



dq 



q 2 n{E q 



x 



/ ( lof 



\m 2 — x(l — x)(E 2 — A; 2 



(32) 



where E a = E p + 0\E g , = |p + q|, the sum is over o = ±, and the angular brackets denote 
the average over the angles of p and q. It is convenient to change the angular integration 
variable to k. The integral over x in the angular average then reduces to 

|m 2 -x(l -x)(E% - k 2 )\\ 



dx ( k>£ 
\ 



1 rp+q 



77T 



2pq J\ 



dkk [f 2 (E a , k) — 2] 



(33) 



p-q\ 



where the function in the integrand is 



h{E,k) 



' E 2 -M 2 \ 1/2 
E 2 -k 2 



log 



(E 2 - k 2 ) 1 ' 2 + (E 2 - Ml) 1 / 2 



' M 2 -E 2 \ 1/2 
E 2 -k 2 



(E 2 — k 2 ) 1 / 2 — (E 2 — Ml) 1 / 2 
at an 



E 2 — k 2 \ 1/2 



Ml- E \ 



k 2 <E 2 -4m 2 , 
E 2 - 4m 2 <k 2 <E 2 



' Ml-E 2 \ 1/2 
k 2 -E 2 



log 



2U/2 



(M 2 k - E 2 ) 1 / 2 + (k 2 -E 
(Ml - E 2 ) 1 / 2 - (k 2 - E 2 ) 1 / 2 



E 2 <k 2 



and Ml = 4m 2 + k 2 . Our final result is 



J(2) 



(47r) 6 U 



- + 2)J( + K 2 



where K 2 is the function of f3m defined by the following integral: 



K, 



?4 I (] l> 



pn(E p ) fP qn(E q) 



E n 



JO Cj n Jp-q 



In the limit m — > 0, K 2 reduces to 

(4vr) 4 



72 



log 



4vrT 



m 



C'(-i) 
C(-i) 



(34) 



(35) 



(36) 



(37) 



where ((z) is the Riemann zeta function. In that same limit, J\ reduces to 

(47T) 2 



Jl 



12 



1 + 2e 



log + 1 + ^-^ 

m C(-l) 



+ 0(t l 



(38) 



Inserting ( |37"1) and (|3~8"D into (^), we obtain the expression for in the limit m —> 0: 



J(2) 



6c 



144(4tt) 2 V4ttT, 

This agrees with the analytic result first obtained by Frenkel, Saa, and Taylor 



(39) 
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E. Three thermal factors 



The integral in fllOD is finite in three spatial dimensions, so we can set e = from 
the beginning. After using the delta functions to integrate over p , q Q , and r , the integral 
reduces to 



(4vr) 6 ./o ^ E p Jo ^ E q Jo E r %\ Ej T -k 2 -m 

where E aT = E p + aE q + rE r , k = |p + q + r|, the sum is over a — ± and r = ±, the angular 
brackets denote the average over the angles of p, q, and r, and V denotes the principal value 
prescription for the poles in the propagator. Before averaging over the angles of p, q, and 
r, it is convenient to use the symmetry in the integration variables to impose the restriction 
r < q < p while multiplying by 3!. We can then average over angles using the identity 

(F(|p + q + r|)> = -L ( P+q+r dkkw(p,q,r,k) F{k) , (41) 
Apqr Jo 

where the weight function for the case r < q < p is 

w(p, q, r, k) = 2k 9(q + r — p) , 0<k<\p — q — r\, 

= k + q + r — p, \p — q — r\ < k < p + r — q , 

= 2r, p + r — q<k<p + q — r, 

= p + q + r — k , p + q — r<k<p + q + r. (42) 

Integrating over k, our final result is 

1(31 = Wf K3Ti ' (43) 

where K 3 is the function of (5m defined by the following integral: 

96 r°° pn(E p ) rp qn(E q ) fi rn(E r ) ^ 
K3 = r4 P ~Ep~ Jo dq —Jo dr ^~^M E ^P + 1 + r) 

-f 3 {E aT ,p + q-r)- f 3 {E aT ,p-q + r) + f 3 (E aT ,p-q-r)} . (44) 
The function in the integrand is 

ME,p) = plog m2 ~ E2 + P 2 + 2 ( m * _ E 2y/2 atan P , E 2 < m 2 , 

m z [mr — E^yi 2, 

_„ 1nr \E 2 -m 2 -p 2 \ 2 2 1/2 (ff 2 -m 2 ) 1/2 +P F 2. m 2 a ^ 
-plog \-{E -m ' lo gT7T^ 2TT^ M E >m . (45) 

The integral in the limit m —>■ was calculated by Frenkel, Saa, and Taylor numeri- 
cally |§ and by Arnold and Zhai analytically ||: 



K, - ^ 



48 



7 + C'(-l) C'(-3) 



15 C(-l) C(-3) 



(46) 



Its numerical value in this limit is K 3 — ^ 453.51 
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1a 2a 2b 

FIG. 1. One-loop, two-loop, and three-loop diagrams contributing to the free energy. 



F. Total 



The final result for the basketball diagram is obtained by inserting (fL^), (|29|), (|35|), and 
1) into ©: 



Ji 



ball 



(47r) 6 Vm 



2 23 35 + tt 2 ' 
+ — + — : + C 



3e 2 



2e 



4 i 

m + 



6 17 
+ 4d 



6 



+ ( - + 12 ) JfT* + (6K 2 + AK 3 )T 



(47) 



To obtain the Laurent expansion including all terms through order e°, it remains only to 
expand the factor (/z/m) 6e and the integral J\ in powers of e. In the limit m — > 0, ( $7\ ) 
reduces to the analytic result obtained by Arnold and Zhai ||: 



X, 



ball 



1 / 



24(4vr) 2 V4vrT 



1 , 91 , C(-1) qC'(-3) 
e 15 + C(-l) C(-3) 



(48) 



III. THREE-LOOP FREE ENERGY 



Using our result for the massive basketball diagram, we can calculate the free energy for 
a massive scalar field theory with a 4 interaction to three-loop order. The Lagrangian for 
the field theory is 

C = \d^4> - \m 2 <f - ^2 W + A£ , (49) 

where AC includes the counterterms. We define the parameters m and g by dimensional 
regularization and minimal subtraction, so they depend implicitly on the renormalization 
scale \x. 



A. One loop 



The free energy at zeroth order in g is given by the one-loop diagram labeled in Fig. |T|: 

-^o = -^X , (50) 
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where the sum-integral 1' Q is defined by (|A.1|) in the Appendix. Keeping only the 
temperature-dependent term, the result for the one-loop contribution to the free energy 
is 



To 



2c 



JoT 



(51) 



2(4tt) 2 \mj 

where Jo is the function of (3m defined in ( |A.5| ). In the limit e — > 0, this function reduces to 

16 '°° 



Jo 



e=0 



3T 4 



poo p* 

/ dp —n(E p 

JO -C/T) 



(52) 



B. Two loops 

The free energy at second order in g comes from the two-loop diagram labeled la in Fig. [1], 
and also from inserting the order-g 2 mass counterterm Aim 2 into the one- loop diagram: 



dT 3 2 



dm 2 



(53) 



The expression for the diagram la is 



(54) 



where the sum-integral X n is defined in (|A.2 ). Keeping only the temperature-dependent 
terms, this diagram is 



F is 



a 



](Att) 2 \m 



Ac 



2 ( - + 1 + 0(e) ) J x m 2 T 2 + J 2 T 4 



(55) 



where at = g 2 /(47r) 2 . The pole proportional to Jifh 2 T 2 is canceled by the last term in (|53|). 
The identity ( |A.6| ) is useful for calculating the derivative dJo/drh 2 in that term. The mass 
counterterm is thereby determined to be 



Aim = — am . 
2e 

Our final result for the two-loop contribution to the free energy is 

a 



-2{L + l)J l m 2 T 2 + J 2 T 4 



(56) 



(57) 



where L = log(/i 2 /m 2 ) and J\ is the function of (3m defined in ( |A.5| ). In the limit e — > 0, it 
reduces to 



Ji 



e=0 



8 f°° p 2 

Y 2 L dp W n{Ep) - 



(58) 
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C. Three loops 



The free energy at second order in g 2 comes from the three-loop diagrams labeled 2a and 
2b in Fig. [l], and also from inserting counterterms into the one-loop and two-loop diagrams: 



T 2 = ^2a + T 2h + +— +-- 



dm 2 



r 



2 (dm 2 )- 



The expressions for the diagrams 2a and 2b are 

1 



2a 



16 



9^2 , 



F2b — — ^ball 



(59) 

(60) 
(61) 



Using the expressions for I n in the Appendix and for Xbaii in fl4"7|), the temperature-dependent 
terms in these diagrams are 



a 



2a 



16(4tt) 2 Vm 



( : k- 



2 ix 2 
e 6 



J 2 m 4 + 2 



1 1 



7T 

~6~ 



2rTi2 



1 



+2 I - + 1 ) J x J 2 m 2 T 2 - -J(T* - J'fJ 2 T 



1 



2rr4 



a 2 



2 b 



16(4tt) 2 Vm 



17 



3e 

+ 4^ J X 2 T 4 



(62) 



(63) 



The identity ( |A.6|) is useful for computing the derivatives with respect to rh 2 in (|5^). The 
pole in e proportional to J 2 T 4 in ( |62"D is canceled by the A^ 2 term in (|59|) . After taking 
into account the terms in ( |59"D involving Axm 2 , the remaining poles in ( |62"D and (p3) are 
proportional to Jifn 2 T 2 and are canceled by the A 2 m 2 term in (|59D . The new counterterms 
that enter at this order are 

a 2 3 



2e 



a# 2 , 



1 



A 2 m = — r am. 

2 V2e 2 24e/ 

Our final result for the three-loop contribution to the free energy is 



(64) 
(65) 



To 



a 



16(4tt) 2 



28 



7T 



4 



(L + l) 2 J 2 m 4 + I 4L 2 + - 4 - — - 1 Jxm 2 T 2 



+2(L + 1) Ji J 2 m 2 T 2 - (3L + 4) JfT 4 - JfJ 2 T* -\2K 2 + -K 3 T 



where C\ = —9.8424. The functions K 2} K 3 , J , and Jj, are given by (p6|) , @), 
(H), and J 2 is 



J 2 



/•oo 1 

4 / dp ^rn(£; p 

JO tin 



(66) 
and 

(67) 



The complete result for the free energy to order o? is the sum of 



and (pq) . 
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1a 2a 2b 

FIG. 2. One- loop and two- loop diagrams that contribute to the self-energy. 



IV. PHYSICAL PARAMETERS 

In this section, we express the free energy in terms of the physical mass m of the boson 
at zero temperature and the physical coupling constant g defined by the threshold scattering 
amplitude at zero temperature. 

A. Physical mass 

The physical mass m of the scalar particle at zero temperature is given by the location 
of the pole in the propagator. If II(P 2 ) is the self-energy function in Euclidean space, then 
m satisfies 

P 2 + m 2 + n(P 2 ) = at P 2 = -m 2 . (68) 

This equation can be solved perturbatively for m 2 as a function of the parameters fh and g 
defined by dimensional regularization and minimal subtraction. To express the free energy 
in terms of m to three-loop order, we need to calculate m 2 to order g 4 . 

The one-loop self-energy IIi, which is independent of P 2 , can be written 

III = IIxa + Axm 2 . (69) 
The expression for the one-loop diagram la in Fig. |2| is 

Ilia = Ifh, (70) 



where the one- loop integral Ji is given in (|A.12 ). Adding the counterterm in (|56|), the 
one-loop self energy is 

III = ~(L + l)arh 2 . (71) 



The two- loop self-energy function, which depends on P 2 , is 

n 2 (P 2 ) = n 2a + n 2b (P 2 ) + H^m 2 + ^A l9 2 + A 2 m 2 . (72) 
The expressions for the two-loop diagrams 2a and 2b in Fig. are 
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(73) 
(74) 



To calculate the physical mass to order a 2 , we need the value of IT 2a (P 2 ) only at P 2 = 
Inserting the values for Ii and / 2 from the Appendix and the value of / sim (-m 2 ) from 
we obtain 
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(75) 
(76) 



Combining all of the terms in (|T2"D, the value of the two- loop self-energy at P 2 



-m is 
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(77) 



The solution to the equation (|68D for m 2 to order a 2 is 



m 



m 2 + 111 + II 2 (-m 2 ) 



Inserting ( [71] ) and ( |7?D into (0), our final result for the physical mass to order a 2 is 
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B. Physical coupling constant 



A convenient physical definition of the coupling constant g is that the amplitude for 
2 — > 2 scattering is exactly — g 2 at threshold where all four particles have four-momentum 
p = (m, 0). To express the free energy in terms of g to three-loop order, we need to calculate 
g 2 to order g 4 . 

The one-loop expression for the negative of the scattering amplitude at threshold is 
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(80) 



where the bubble integral is defined in (|3Cf ). Using the result (pl|) , the values of the bubble 
integrals that appear in (|80|) are 
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(81) 
(82) 



We have neglected the difference between m and to in the bubble integrals because it is 
higher order in a. Inserting 



81Q and (|82|) together with ( p4|) into (pOf) , our final result for 



the physical coupling constant a = g 2 /(47r) 2 is 
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(83) 



C. Three-loop free energy 



To express the three- loop free energy in terms of the physical mass and coupling constant, 

insert them 



we need to invert (|79|) and (jS3|) to obtain m 2 and g 2 in terms of m 2 and g 1 
into our expression for the free energy, and expand to order a 2 . Inverting ([FjJ 
obtain 



and 
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where L = log(/i 2 /m 2 ), not to be confused with L = log(/i 2 /m 2 ) in (|79| ) and fl8"3|). Upon 



inserting these expressions into the sum of (|51~D , (|57|) , and (|66l) , and expanding to order ft 2 , 
all the terms that depend on L cancel. Our final expression for the temperature-dependent 
contribution to the free energy in terms of physical parameters is 



1 



2(4tt) 
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2J 2 + J 2 J 2 + 2K 2 + -K 3 



(86) 



where K n and J n are the functions of (3m defined by (|36|) , (|J), d52|) , (pq) , and (p7|) . This 
expression is remarkably compact. 

The effect of the interaction on the free energy (which is the negative of the pressure) 
is illustrated in Fig. |^. We normalize the free energy to that of an ideal gas of particles 
with the same physical mass m, which is given by T§ in ([Sip. We plot J r /J r o as a function 
of T/m on a log scale for two different values of the physical coupling constant: a = 0.1 
and ft = 0.4, which correspond to g = 3.97 and g = 7.95, respectively. The dashed lines 
are the free energies truncated after the order-ft terms. The solid lines are the free energies 
truncated after the order-ft 2 terms. 

For T <C to, the three- loop result for the free energy (|86|) approaches 
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(87) 



The exponential approach to the free energy of an ideal gas is evident in Fig. |3|. Note that 
the order-ft 2 correction is smaller than the order-ft correction only if ft < |. For T ^> to, 
(BGI) approaches 
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FIG. 3. Free energy normalized to that of an ideal gas vs. T/m for a = 0.1 and 0.4. The 
dashed and solid lines are the two- and three-loop approximations to the free energy, respectively. 

In the order-a 2 correction, the linearly divergent T/m term is the first of a series of infrared 
divergent terms that behave like a n+1 (T / m) 2n . These terms come from the ring diagrams 
which, when summed to all orders, give a correction of + (5v^6/3)a 3//2 . The logarithm in 
the order-a 2 correction term arises from the running of the coupling constant. It can be 
absorbed into the order-a correction term by replacing the physical coupling constant a 
by a(T), the MS coupling constant with renormalization scale \i — T . For T ^> m, we 
expect the three-loop result to be a good approximation only if the a 2 T/m correction is 
small compared to the a correction, which requires T <C m/{na). 



V. SUMMARY 



We have reduced the thermal basketball diagram for a massive scalar field theory with 
a 4 interaction to three-dimensional integrals that can be evaluated numerically. As an 
application, we calculated the free energy for this theory to order a 2 . The result is par- 
ticularly simple if the free energy is expressed in terms of the physical mass and coupling 
constant. Another useful application of our result for the massive thermal basketball dia- 
gram would be to extend the calculation of the free energy for the massless theory using 



screened perturbation theory to three-loop accuracy [16 
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APPENDIX A: ONE-LOOP SUM-INTEGRALS 



The one-loop sum-integrals required to calculate the free energy to order g 4 are 

?o = ~f p log (P 2 + m 2 ) , (A.l) 

(A.2) 



T 

- L - ,r n 



p (P 2 + m 2 ) n 



The sum-integral X' is the derivative of T n with respect to its index evaluated at n — 0. 
These integrals satisfy 
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dm? 
d 



j n 

dm 2 n 



-nX, 



■n+l 



(A.3) 
(A.4) 



The specific sum-integrals that are required are I' , I±, and I 2 . The temperature- 
dependent terms in the sum-integrals can be conveniently expressed in terms of the following 
integrals: 



4eT(!) 
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(A.5) 



These integrals are functions of f3m only and satisfy the recursion relation 

d 

m—J n = 2eJ n - 2(f3m) 2 J n+1 . 
dm 



(A.6) 



If we separate out the temperature-dependent terms in the one-loop sum-integrals, the 
resulting expressions are 
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(A.7) 
(A.8) 
(A.9) 
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To calculate the physical mass and coupling constant, we also need the one- loop Eu- 
clidean momentum integrals I\ and I2 defined by 



p (P 2 + m 



2\n 



(A.10) 



These integrals satisfy 







dm 2 



-nl, 
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(A.ll) 



The integrals Ii and I2 are identical to the temperature-independent terms in ( |A.8| ) and 
respectively. Expanding around e = 0, these integrals are 
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Addendum 



After this paper was completed, J.-M. Chung (DJ provided us with analytic expressions 
for the coefficients Co and C\ defined in (|1^) and (|28|): 

+ Y C(2)-2C(3), 



Ci = - 



59 



C(2). 



These analytic expressions can be derived using the methods described in Refs. [ H| H^ ] . 
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